Reduced Visibility of Rabi Oscillations in Superconducting Qubits 
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Coherent Rabi oscillations between quantum states of superconducting micro-circuits have been 
observed in a number of experiments, albeit with a visibility which is typically much smaller than 
unity. Here, we show that the coherent coupling to background charge fluctuators [R. W. Sim- 
monds et al, Phys. Rev. Lett. 93, 077003 (2004)] leads to a significantly reduced visibility if the 
Rabi frequency is comparable to the coupling energy of micro-circuit and fluctuator. For larger 
Rabi frequencies, transitions to the second excited state of the superconducting micro-circuit be- 
come dominant in suppressing the Rabi oscillation visibility. We also calculate the probability for 
Bogoliubov quasi-particle excitations in typical Rabi oscillation experiments. 
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PACS numbers: 85.25.Cp,03.65.Yz,03.67.Lx 

Possible applications in quantum information process- 
ing have renewed the interest in quantum coherent phe- 
nomena of micrometer-scale Josephson junction (JJ) cir- 
cuits (Ij. During the past years, several experiments have 
demonstrated Rabi oscillations of a macroscopic vari- 
able including chargeH, phase a combinations of 
both 4,1 M I ^ii'i flux [3 , which persist up to several mi- 
croseconds 0. However, in many of the reported ex- 
periments, the Rabi oscillation visibility is significantly 
smaller than unity even at times short compared to the 
decoherence time. Characteristic values are as small as 
10% (Ref. ^ and 50% (Refs. H 0), respectively, which 
indicates either substantial leakage out of the computa- 
tional basis or a pronounced randomization of the state 
occupation over short time-scales. In the following, we 
discuss three possible mechanisms that result in a re- 
duced visibility, namely (i) coherent coupling to a back- 
ground fluctuator which induces transitions between the 
eigenstates of the superconducting micro-circuit 7J; (ii) 
population of higher excited states caused by a strong 
driving field, i.e., the AC current or voltage applied to 
drive Rabi oscillations; and (iii) the excitation of Bogoli- 
ubov quasi-particles. 

Background charge fluctuators are presently consid- 
ered one of the dominant source of decoherence in su- 
perconducting qubits (squbits) [ElMHiOl- However, 
the theoretical analysis of fluctuators has so far been fo- 
cused on fluctuators with an incoherent dynamics which 
are coupled to the squbit by an Ising-like interaction. Mo- 
tivated by the results of Refs. 0and[ll we show that the 
coherent dynamics of a squbit-fluctuator system with a 
transverse exchange coupling leads to a reduced visibility 
and beating in the Rabi oscillation signal of the squbit. 
For Rabi oscillation frequencies large compared to the 
squbit-fluctuator coupling strength J, the fluctuator does 
not provide an efficient mechanism for visibility reduc- 
tion on short time-scales. We show that, in this regime, 
coupling to the second excited squbit state suppresses the 



visibility to 0.7 for Rabi frequencies larger than 150 MHz 
in the experiment in Ref. 13- This mechanism persists 
for adiabatic switching, in stark contrast to the leakage 
scenarios discussed previously 0, ■ Finally, we cal- 
culate the excitation rate of Bogoliubov quasi-particles 
and show that the quasi-particle excitation probability is 
small as long as the JJ is in a zero-voltage state. 

Coherent squbit-fluctuator coupling. - While fluctu- 
ators, trapped charges in bi-stable potentials, have so 
far mainly been considered as source of decoherence for 
squbits, recent experiments provide strong evidence that 
some fluctuators are coupled coherently to a squbit. An 
important conclusion drawn from the anti-crossings in 
the microwave spectra of a phase-squbit and the sub- 
sequent observation of real-time oscillations is that 
the squbit-fluctuator coupling has a transverse compo- 
nent rather than the Ising-form considered in previous 
theoretical analysis 0, llll |. Because of the transverse 
coupling component, fluctuators induce transitions be- 
tween different squbit basis states, randomizing the oc- 
cupation probabilities. In order to analyze the dynam- 
ics of a coupled squbit-fluctuator system, we introduce a 
pseudo-spin notation and deflne ~ (|1)(1| — |0)(0|)/2 
and Sx = (|1)(0| + |0)(l|)/2 in terms of the ground (flrst 
excited) state |0) (|1)) of the squbit, such that js^ =J,) 
represents the squbit ground state. Similarly, we de- 
flne the pseudo-spin states \Iz =|) and \Iz =T) as the 
ground and first excited states of the fluctuator, respec- 
tively. Rabi oscillations between the squbit basis states 
are driven by an AC current in resonance with the level 
splitting huiQ of the squbit. In our pseudo-spin notation, 
the AC current acts as transverse magnetic field and, 
in the co-rotating frame. 
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where S = h{ujeg — ujio) is the detuning of the fluc- 
tuator level splitting hweg relative to the squbit. The 
ansatz for the transverse squbit-fluctuator exchange cou- 
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pling with couphng strength J is microscopically moti- 
vated from the obseryation that the critical current of 
a JJ may depend on the position of a fluctuator in its 
bi-stable potential 0. As shown below, a fluctuator 
influences the squbit dynamics strongly ii \S\ < J or 
|(5 ± 6^1 < J. For the J J in Ref. J/h ~ 25 Mhz 
while the typical level spacing between different fluc- 
tuator resonances is of order 60 MHz. This allows us 
to restrict our analysis to one fluctuator with minimum 
\S\ or |5 ± 63^1 first. For the experimental temperature 
T ~ 20 mK <C huJiQ/kB,fiujeg/kB, at the beginning of 
the Rabi pulse, t = 0, both squbit and fluctuator are 
in their ground-state and IV'(O)) = | i; i) in the product 
basis \sz] Iz) of the squbit-fluctuator system. The exper- 
imentally accessible quantity is the probability pi{t) for 
the squbit to occupy its first excited state as a function 
of Rabi pulse duration, pi{t) = X]/^=t j. Kt; 
While the energy level splitting hujeg of a fluctuator is 
fixed, the squbit level splitting Hwiq can be tuned via 
the DC bias current through the J J, which allows one to 
measure squbit Rabi oscillations for varying 5 and hx- 

We now calculate the dynamics of IV'(O)) = | i; i) and 
Pi{t) as a function of 5 and bx- The squbit-fluctuator 
exchange coupling gives rise to a linear dependence of the 
eigenenergies on J for |5| < J and |5 ± < J, where 
cross-relaxation processes between squbit and fluctuator 
reduce the visibility of the squbit Rabi oscillations 
While pi{t) is readily obtained from integration of the 
Schrodinger equation for arbitrary bx, S, and J, here we 
focus on the two cases 5 = and 5 ± bx — 0, where the 
cross-relaxation between squbit and fluctuator is most 
efficient. For 6 = 0, we find 



(5=0 



-[f-cos(-)cos( V ) (2) 

sm(f )sm( ^ 2;^ ) i 
V^l + {2bx/JY ^ 

l[l_cos(|)cos(^)]. 



For 5 ±bx — Q, the exact expression for p\{t) is too long 
to be presented here, but for \bx\/ J > 2 it is well approx- 
imated by 



Pi{t)\ 



Equations (01 and lj3l describe the Rabi oscillations of 
a squbit in the presence of a fluctuator. The trans- 
verse coupling to a fluctuator introduces an additional 
Fourier component which leads to beating of the Rabi- 
oscillation signal, in agreement with experimental results, 
where typically \bx\/J > 2. For short time-scales, the 
J-dependent factor in Eqs. (O and leads to a de- 
crease in Rabi oscillation visibility. In particular, the 
first maximum in pi [t) is reduced relative to unity. For 



\bx\/ J > 2, we obtain pi(i = irh/bx) — 1 — {irJ/Abx)'^ and 
pi{t = T^fi/bx)^ 1 - (TTJ/Sbx)^ tor 6 = and S±bx = 0, 
respectively While a resonant fluctuator reduces 

Pi{t = Tvh/bx) to 0.85 for \bx\/J — 2, Rabi oscillations 
with amplitude I are predicted to emerge for \bx\/ J ^ f , 
in stark contrast to experimental results where the first 
local maximum of pi (i) is of order 0.5 even for large bx ■ 
In a realistic system, decoherence of the squbit and 
fluctuator dynamics will damp out both the Rabi oscilla- 
tions and the beating predicted for pi{t) in Eqs. lO and 
In order to discuss reduced visibility in the pres- 
ence of decoherence, we focus on a simple model system 
in which the fluctuator has a finite decoherence rate 7//1 
while its relaxation rate vanishes ('diagonal coupling' of 
a bath and fluctuator). The dynamics of the system is 
determined by the Master equation for the density ma- 
trix [13, 



(4) 



A similar system was analyzed in Ref. 18 using different 
numerical and approximate analytical approaches. Here, 
we provide an analytical solution of Eq. for \bx\/ J ^ 1 
where, to leading order in J, _ff ~ 51^ + Jsxix + bxSx- 
Then, the sixteen coupled differential equations in Eq. Q 
decouple into sets of four differential equations. We cal- 
culate pi{t) = ^1 j^(t; /z|/5(i)| VjIz) from the explicit 
solution of Eq. Q. For J = and p(0) | i; i |, 
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+ (1 + ^^_)e-(7-V^^)*/2n]g-^b.tM\ 

Note the intriguing dependence of p\{t) on the fluctu- 
ator decoherence rate 7//1 and the coupling strength J 
(Fig.HJ). For 7/J< 1, 

PiW-^[l-e-'^*/^''cos(g)cos(^)] (6) 

shows the beating already derived in Eq. |5J|, but with a 
finite damping rate ^/Ih. In contrast, for 7/ J ^ 1, 



g-4J=tA7eos(Y)] 



(7) 



exhibits single-frequency oscillations with a damping rate 
4J^/ft7 (see also R,ef. Il8r. This can be understood from 
the obseryation that, for 7 ^ J, the phase of the fluctua- 
tor is randomized on a time-scale short compared to h/ J 
and the mean field J{Ix) acting on the squbit is averaged 
out to leading order, which restores the single-frequency 
oscillations in Eq. jT)). 

We discuss next to which extent far off-resonant fluc- 
tuators with \S\ » l^^;! » J reduce the Rabi oscil- 
lation visibility. Because of the perturbative parame- 
ter J/\S\ ^ 1, the visibility reduction in pi{t) ~ (1 — 
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FIG. 1: (color online). Rabi oscillations for a squbit- 
fluctuator system. The probability pi{t) to find the squbit 
in state |1) is obtained from numerical integration of Eq. @ 
(solid line) and the analytical solution Eq. Q (dashed line) 
which is valid for bx/J S> 1. For weak decoherence of the 
fluctuator, 7/J < 1 [(a) and (c)], Pi{t) shows damped beat- 
ing. For 7/J > 1, damped single- frequency oscillations are 
restored. The fluctuator leads to a reduction of the first max- 
imum in pi{t) to ~ 0.8 [(a) and (b)] and ~ 0.9 [(c) and (d)], 
respectively. The parameters are (a) bx/J ~ 1.5, 7/J = 0.5; 
(b) bx/J = 1.5, 7/J = 1.5; (c) bx/J = 3, 7/J = 0.5; (d) 
&./J = 3, 7/J = 1.5. 



J^/4(5^) sm'^{bxt/2h) is small for a single ofF-resonant fluc- 
tuator. More generally, a large number of ofF-resonant 
fluctuators interacts with the squbit, and the problem is 
closely related to the dynamics of an electron spin cou- 
pled to a bath of nuclear spins by the hyperfine contact 
interaction To leading order in the coupling con- 

stants, a set {i} of off-resonant fluctuators with coupling 
constants Ji and detunings 6i reduces the Rabi oscilla- 
tion visibility by Jf/ASf. Introducing the distribution 
hmction P{E) of the fluctuator level splittings, for con- 
stant J, = J, 



E 



E-hu:iQ\>\b^\;E>E, 



dEP{E) 



(8) 



where the integral is evaluated for all energies with 
\E — /lo^iol > 1 I (off- resonance condition) larger than 
a lower energy cut-off E(, given by the fluctuator de- 
coherence rate. Evaluating the integral for the distri- 
bution function P{E) cx \/E characteristic for fluctu- 
ator level spacings (l//-noise), we find ^iJ'^/'^5f ~ 
2J'^P{hjidiQ)/A\bx\ — 0.2J/|6a;| for the experimental pa- 
rameters of Ref. [7|. If cross correlations between fluctu- 
ators can be neglected, the reduction of the first maxi- 
mum in pi (t) is obtained by summing the contributions 
from resonant fluctuators with i5 ~ and (5 ± &a; ~ 
[Eqs. lEJ and ©] and off-resonant fluctuators [Eq. ©]. 
An upper bound is given by 1 — pi(i = irh/bx) < 
2J'^P{hujio)/4\bx\ + i3/32){TrJ/bx)'^, which decreases with 



increasing \bx\- This shows that, even when both reso- 
nant and off-resonant fluctuators are taken into account, 
the small visibility pi{t — iih/bx) < 0.5 observed in the 
limit of large Rabi frequencies, [fc^l/Zi ^ J/h in Ref. 
cannot be effected by fluctuators. 

Energy shifts induced by AC driving field. - Exploring 
the visibility reduction at a time-scale of 10 ns requires 
Rabi frequencies \bx\/h > 100 MHz. We show next that, 
in this regime, transitions to the second excited squbit 
state lead to an oscillatory behavior in pi (t) with a visi- 
bility smaller than 0.7. For characteristic parameters of 
a phase-squbit, the second excited state 1 2) is energeti- 
cally separated from |1) by 1^21 = 0.97^10 u^- Similarly 
to |0) and |1), the state |2) is localized around the local 
energy minimum in Fig. |2Ia). For adiabatic switching 
of the AC current, transitions to |2) can be neglected as 
long as \bx\ ^ hAu = ^{ujio — u)2i) — O.OSftwio. How- 
ever, for bx comparable to hAu, the applied AC cur- 
rent strongly couples |1) and |2) because (2I0I1) 7^ 0, 
where (fi is the phase operator. For typical parameters, 
bx/hAuj ranges from 0.05 to 1, depending on the irradi- 
ated power [3, 01 . Taking into account the second ex- 
cited state of the phase-squbit, the squbit Hamiltonian 
in the rotating frame is 14] 

H = -hAuj\2){2\+bx (|0)(1| + V2\l){2\ +H.c) /2. (9) 



In the following, we neglect fluctuators, which is valid for 
the short-time dynamics if [fo^;! S> J. 

The time-evolution of ['(/'(O)) — |0) is readily calculated 
by integration of the Schrodinger equation. Expanding 
IV'(t)) in {bx/2hAuj), we find that 



_ 3 / bx y 

2\2hALo) 



n 2 



sm 



bx [l-{bx/2hALof]t 



(10) 

exhibits single-frequency oscillations with reduced visi- 
bility [Fig. |2Ib)] . Part of the visibihty reduction can be 
traced back to leakage into state |2). More subtly, off- 
resonant transitions from |1) to |2) induced by the driving 
field lead to an energy shift of |1), such that the transition 
from |0) to |1) is no longer resonant with the driving field, 
which also reduces the visibility. For bxl2hAu) = 1/3, 
corresponding to bx/h = 150 MHz in Ref. 0, Eq. ifTUI) 
predicts a visibility of 0.7. Note that the off-resonant 
transition to |2) discussed here is induced by a large am- 
plitude of the driving field and not by additional frequency 
components that result, e.g., from non-adiabatic switch- 
ing of the driving field | l3| or current noise in the micro- 
circuit H^. Rather, the visibility reduction in Eq. ifTUI) 
corresponds to a steady-state solution of the Schrodinger 
equation and persists even for adiabatic switching [2l|. 
Whether such off-resonant transitions are responsible for 
the reduced Rabi oscillation visibility can be tested ex- 
perimentally from the dependence of the Rabi frequency 
and amplitude on bx- Equation H10|l predicts a non- linear 
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FIG. 2; (color online), (a) Level scheme for the squbit in 
Refs. 0, lifl . (b) pi (t) obtained from numerical integration of 
the Schrodinger equation (solid line) in comparison with the 
analytical result in Eq. IIUII (dashed line) for bx/hAuj = 1/3. 



dependence of the Rabi frequency on , while the visi- 
bility is predicted to decrease quadratically with ■ De- 
creasing visibility with increasing driving field has indeed 
been observed in some Rabi oscillation experiments |22l| . 
but further quantitative analysis is required to determine 
the explicit functional dependence. 

Excitation of quasi-particles.- Transient high- 
frequency components in the switching pulses induce 
leakage to states outside the computational basis. Bo- 
goliubov quasi-particle excitations represent a large class 
of excited states which is often ignored for the discussion 
of leakage because the excitation gap is large. For both 
charge-based and phase-based squbits, the excitation 
of quasi-particles during the measurement process has 
been quantified and is known to limit the decoherence 
time j^lUISl- I'^ phase-squbits, quasi-particles trigger 
leakage to finite-voltage states or to the excited squbit- 
states 1 2) or |3) by incoherent tunneling across the J J at 
a rate determined by the single-particle tunneling rate 
and the total number of quasi-particles Q . 

While the number of quasi-particles created during 
read-out can be decreased by increasing experimental 
waiting times j2^, Bogoliubov quasi-particles are also 
excited during an experimental cycle by high-frequency 
Fourier components of an external current or voltage 
pulse with u} > 2A/h, where A is the energy gap of 
the superconductor. Microscopic mechanisms for quasi- 
particle excitation include, e.g., dissipative tunneling 
through the JJ. We calculate the excitation rate F of 
quasi-particle pairs for an oscillating current component 
6Iu}C0s{ujt) with LO > 2A/H. Using a semiclassical ap- 
proximation, the high-frequency current induces oscilla- 
tions S(j){t) of the superconducting phase around its equi- 
librium value (jjQ — arcsin(///c) ~ ti'/Z, where / and 
Ic are the DC bias current and critical current of the 
JJ, respectively. Amplitude and frequency of S(j){t) are 
determined by the externally imposed current, S(j){t) = 
— [{eSIuj)/{huj'^C)] cosojt, where C is the JJ capacitance. 
Dissipative tunneling through the J J - processes in which 



a quasi-particle pair is excited under absorption of en- 
ergy huj from the bias current - is described by both the 
normal current and the Josephson cosine-term j2^. To 
lowest non- vanishing order in (eSI,,,) / (huj'^C), the time- 
averaged power dissipation P at the JJ is |25( 



P = 



1 fSI^ 
2\Cu; 



huj 



[In {huj) + hihuj) COS ( 



1 /SI^V e 



(11) 



with the standard expressions for the normal current /„ 
and the Josephson cosine-term l2- For huj — > 2A+, /„ 
is determined by the normal-state resistance i?„ of the 
J J |2fii], /n(2A+) = (7r/4i?„)(2A/e). The second line of 
Eq. Hll|l is valid for squbits biased close to the critical 
current, where 0o — 7r/2. The quasi-particle excitation 
rate is obtained from P via |27| 



P 



(12) 



F is negligibly small unless C < 1 fF and 61^; > 1 pA. 
With typical parameters i?„ = 29 fl, C — 1 fF, and 
Lj — AKks/h, we obtain F = Imsec"^ x ((5/c^/pA)^. 
For instantaneous switching, a current pulse with am- 
plitude lac ~ 10 nA and duration T = 10 usee, 5Iuj — 
lac/^T = IpA for Lo = AKks/fi- These values show that 
quasi-particle excitation by dissipative tunneling cannot 
explain the substantial visibility reduction evidenced in 
current Rabi oscillation experiments for phase squbits. 
However, the excitation of Bogoliubov quasi-particles is 
relevant in view of the ultimate goal of reducing quantum 
gate errors to less than 10~^. 
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